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Abstract 



We study theoretically the ultrafast nonlinear optical response of quantum 
well excitons in a perpendicular magnetic field. We show that for magnetoex- 
citons confined to the lowest Landau levels, the third-order four- wave-mixing 
(FWM) polarization is dominated by the exciton-exciton interaction effects. 
For repulsive interactions, we identify two regimes in the time-evolution of the 
optical polarization characterized by exponential and power law decay of the 
FWM signal. We describe these regimes by deriving an analytical solution 
for the memory kernel of the two-exciton wave-function in strong magnetic 
field. For strong exciton-exciton interactions, the decay of the FWM signal 
is governed by an antibound resonance with an interaction-dependent de- 
cay rate. For weak interactions, the continuum of exciton-exciton scattering 
states leads to a long tail of the time-integrated FWM signal for negative time 
delays, which is described by the product of a power law and a logarithmic 
factor. By combining this analytic solution with numerical calculations, we 
study the crossover between the exponential and non-exponential regimes as a 
function of magnetic field. For attractive exciton-exciton interaction, we show 
that the time-evolution of the FWM signal is dominated by the biexcitonic 
eff'ects. 
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I. INTRODUCTION 



During recent years, the role of many-body correlations has been a dominant theme in 
the studies of the ultrafast nonlinear optical response of semiconductors in the coherent 
regime.0 In the lowest order in the optical fields, the nonlinear optical polarization comes 
from the excitation of two electron- hole (e-h) pair states. The Coulomb interactions between 
these elementary excitations govern the spectral and temporal properties of the nonlinear 
optical response.ia 

In noninteracting systems, such as atom-like two-level systems, the only source of opti- 
cal nonlinearity is the Pauli blocking. This introduces local correlations, due to the phase 
space filling, between the carriers excited by the pump and the probe optical pulses, which 
are separated by a time delay r.@ The Coulomb interaction brings about new nonlinear- 
ities by changing the dynamics of the photoexcited e-h pairs. The simplest theoretical 
approach that accounts for the interactions between photoexcited e-h pairs is the time- 
dependent Hartree-Fock approximation, which leads to the well known semiconductor Bloch 
equations (SBE).li^li The latter have been successfully used to describe a variety of ex- 
perimental results, including, e.g., the ac Stark effect and Rabi oscillations,l'0ll the fast 
exponential decay of the time-integrated (TI) four-wave-mixing (FWM) signal for negative 
time delays (r < O),0lli the time-dependence of the time- resolved (TR) FWM signa]0Hli, 
the photon echoes coming from the continuum of states, and some aspects of phase 
measurements.Ellill Being a mean-field approach, the SBE's replace the Coulomb interac- 
tions by a spatially-averaged effective field of all the photoexcited excitons, and thus neglect 
correlations due to exciton-exciton scattering and bound biexciton states. 

More recent experiments indicated the importance of Coulomb correlations, e.g., in the 
dependence of the nonlinear response on the polarization of the optical fields. Manyof the 
experimental results have been explained in terms of excitation-induced dephasing,EjE2l and 
exciton-exciton scattering processes.ii^ The importance of bound biexciton states in FWM 
was also pointed out, e.g., in Refs. |26|-p2[ Exciton-exciton correlations have been shown to 
play a sig niflcant role in the ac StaVkSecti and in qnantnm beats.B Several systematic 
methods for incorporating many-particle correlations in the coherent nonlinear optical re- 
sponse have been developed,^ including the effective exciton Hamiltonian approach,!!^'! 
the coherently controlled truncation scheme^S^S, the correlation function approach,&@ 
and the canonical transformation approach.ElEl 

In this paper, we study the role of exciton-exciton correlations on the ultrafast nonlinear 
optical response of quantum well (QW) excitons in a perpendicular magnetic field B. The 
effect of a magnetic field on the photoexcited carriers is twofold. First, the conduction 
and valence bands collapse into discrete sets of highly degenerate Landau levels (LL). The 
magnetic field squeezes the electron and hole wave functions (magnetic confinement). The 
relative importance of the Coulomb versus magnetic confinement in determining the exciton 
states depends on the ratio Uc/2Eb = (as//)^, where and Eb are the exciton Bohr 



radius and binding energy at zero field, while / = \Jl/eB and Uc are the magnetic length 
and the cyclotron frequency, respectively (we set ^ = 1). With increasing magnetic field, the 
crossover from Coulomb-bound excitons to magnetoexcitons occurs for I ^ as, where the 
binding energy of the 2D magnetoexciton, Eq ~ e^//, exceeds the Bohr energy Eb ~ e^ /qb- 
In typical GaAs/AlGaAs quantum wells, this crossover to the strong magnetic field regime 
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occurs for S > 4 T. 

The second effect of the magnetic field is that, by freezing the kinetic energy of the pho- 
toexcited carriers in two spatial dimensions, it increases the relative strength of the exciton- 
exciton interactions, i-i In bulk semiconductors, FWM experiments showed a strong en- 
hancement of the FWM signal with increasing magnetic field.0^ At the same time, the 
time-integrated signal exhibited a long non-Markovian decay for negative time delays.ElrH 

In QW's, the exciton-exciton interactions are expected to dominate the coherent optical 
dynamics as the magnetic field increases. When the LL separation exceeds the magnetoex- 
citon binding energy, uJc/\Eq\ ~ a^// ^ 1, the inter-LL transitions are suppressed and the 
optical response is determined by magnetoexcitons confined to the lowest LL (LLL). Note 
that, in an ideal 2D system, the interactions between magnetoexcitons with zero center of 
mass momentum are suppressed due to the electron-hole symmetry.i§H£§ In realistic QW's, 
however, this symmetry is lifted due to, e.g., the differing band offsets, lateral confinement, 
and disorder .ii As was first pointed out in Refs. the e-h asymmetry plays an important 

role in the nonlinear response of magnetoexcitons in QW's. 

We investigate here the role of four-particle correlations in the coherent ultrafast dy- 
namics of 2D magnetoexcitons confined to the LLL. Using the canonical transformation 
formalism0@ii, we show that, due to the discrete nature of LL's, the third-order FWM po- 
larization is determined by a momentum-independent two-exciton amplitude x{^)i satisfying 
a simple equation 

ida{t) - 2Qox{t) = V^^V(t) + F{t), (1) 

where Qq is the pump detuning from the magnetoexciton level. The first source term in 
the rhs represent the two pump- induced polarizations, p{t), interacting via the Hartree-Fock 
potential V^^, while F{t) originates from the exciton-exciton correlations. 

We have calculated F{t) numerically and showed that, for a wide range of parameters, 
the results can be described in terms of an analytical solution derived for strong magnetic 
field. In particular, x(t) can be expressed in terms of the Hartree-Fock solution x^^{t) 
[obtained by setting F(t) = in Eq. (0)] as follows: 

/CO 
dt'S{t-t')x''^{t'), (2) 
-oo 

where S{t) represents a memory kernel whose time dependence is governed by the exciton- 
exciton correlations. For sufficiently strong magnetic fields, such that the magnetic length I, 
is smaller than the chacteristic range of the exciton-exciton potential, a, we have derived a 
simple analytic expression for the memory function S{uj) in the frequency domain, that in- 
corporates both the biexcitonic and the exciton-exciton scattering effects nonperturbatively. 

For attractive exciton-exciton interaction, the biexciton effects dominate. For / -C a, 
when the momentum exchange between excitons is suppressed, we recover the average po- 
larization model (APM). EIH For weaker fields, we show that the time-evolution of the 
FWM signal is determined by both biexcitonic and exciton-exciton scattering effects. For 
short pump durations, we reproduce the biexciton oscillations in the TI-FWM signal. § As 
the pump duration becomes longer than the biexciton binding energy, we find a crossover 
between two exponential regimes in the SR- and TI-FWM spectra, originating from different 
dephasing rates of the exciton and biexciton bound states. 
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For repulsive exciton-exciton interactions, we find two regimes in the time dependence 
of tlie FWM signal governed by a dimensionless parameter G ~ Uga'^ / EqI"^ , where £"0 is 
the magnetoexciton binding energy and Uq is the characteristic strength of the exciton- 
exciton interactions {uq < for repulsive interactions). For |G| 3> 1, the memory function 
is dominated by an antibound resonance, made up from the continuum of the exciton- 
exciton states, with the width ~ uoGe"'*^'. In this regime, the TI-FWM signal exhibits 
oscillations with a period determined by the energy of an antibound state, together with an 
overall exponential decay characterized by the interaction-induced time 7^^. For |G| <C 1, 
the long-time behavior of the FWM signal is determined by the singular (logarithmic) low- 
frequency behavior of S{uj)] this leads to a power low decay of the TI-FWM signal for 
negative time delays: TI-FWMoc (up to log factor). In the intermediate case \G\ ~ 1, the 
crossover between the two regime as a function of time delay occurs on the scale |r| ~ |mo|~^- 
By combining our analytic solution with numerical calculations, we study the crossover 
between exponential and non-exponential regimes as a function of magnetic field. Our 
results demonstrate that four-particle Coulomb correlations lead to a qualitatively different 
dynamics of the two-exciton excitations as compared to the one-exciton states. 

The paper is organized as follows. In Section 0, we summarize relevant facts regarding 
2D excitons in a strong magnetic field. In Section |ITT| we present an expression for the 
FWM polarization obtained in Appendices |^ and |B[ In Section 0, we derive an analytical 
solution for the two-exciton amplitude in the case of strong magnetic field and in Section |^ 
we analyze the memory function in all the regimes of interest. In Section [V^ , we presents 
the results of numerical calculations for the FWM polarization and compare them with our 
analytical solution. Section |VII| concludes the paper. 



II. 2D EXCITONS IN STRONG MAGNETIC FIELD 

In this section, we summarize some basic facts regarding 2D magnetoexcitons.&^i We 
consider the case when the magnetic field is sufficiently strong and the central optical fre- 
quency is tuned to transitions between the lowest Landau levels of the conduction and 
valence bands. In this case, the inter-LL transitions are suppressed so that the photoexcited 
e-h pairs are restricted to the LLL. In the Landau gauge, A = {0,Bx), the electron and 
hole wave-functions, ipk and tpk, are characterized by the y-component of the momentum, 
k = —Xq/I"^, where Xq is the x-coordinate of the cyclotron orbit center, 

Mr) = e^^^-(^-^'='^)^/2^ Mr) = C.(r), (3) 

{L is the system size). The normalized one-exciton wave function, with total momentum p, 
can be expressed in terms of the electron and hole wave functions asil^ 

%ir,,r,) = iV-^/^Ee-^'^^'Vp,/2+.(ri)V5p,/2-.(r2) = ''e^r>■n-^Xy/l^-ir^l^^..)V.l^^ (4) 

k ^ 

where r = ri — r2, R = (ri + r2)/2 are the relative and the average coordinates, respectively, 
= L^/27rZ^ is the LL degeneracy, and z is the unit vector in the direction of the magnetic 
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field. The magnetoexciton creation operator is given in terms of tlie electron and hole 
operators aj, and bl, by the similar expression, 

k 

Consider first the spectrum of a single magnetoexciton. In the one-exciton basis, the 
Hamiltonian H^^^ that describes the interaction between an electron and a hole, is diagonal: 

^S^ = E^p4^p' (6) 

p 

where Ep is the 2D magnetoexciton energy, given by 

'^ = -1 (0,^,e-^'^/^-'('^xp)-^ = -^/f e-^'^/^/o(p^/V4). (7) 

Here Vg = 27re^/fi;g is the 2D Coulomb potential (k is the dielectric constant) and Io{x) is the 
Bessel function. For small center of mass momenta, pi <^ 1, the magnetoexciton spectrum 
is parabolic, 

2 / — 

where Eq = ^wf is the magnetoexciton binding energy and M is the magnetoexciton 

effective mass. For large momenta, pi 3> 1, Eq. (|^) yields Ep = —e'^/upP. Since for strong 
magnetic field the exciton momentum is proportional to the e-h separation, p ~ r/l"^, for 
large momenta we recover the Coulomb potential between a well-separated electron and hole, 
—e'^/Kr. On the other hand, for small momenta, the electron and hole sit on top of each 
other forming a neutral magnetoexciton with the effective mass M.&^l The one-exciton 
Hamiltonian is Hx = Hq + H^,]^ , where Hq describes a noninteracting e-h pair with energy 
Eq + Eq + Eg — uo = Eg — uq, where Eq and Eq are the electron and hole LLL energies, 
respectively. Eg is the semiconductor bandgap, and luq is the pump central frequency (we are 
working in the rotating frame). The eigenstates of Hx, |p) = c?p|0), satisfy the eigenvalue 
equation 

Hx\-p) = ^p\p), Vtp = Eg + Ep- Uq. (9) 

Note that the linear response in strong field is similar to that of a two-level system. This 
can be seen, e.g., by noticing that the optical transition operator, is proportional to the 
exciton creation operator [setting p=0 in Eq. (j^)], 

U^^Y.^ib'-k = N"'d,, (10) 

k 

SO that HU'^Y^) = N^^'^UxdllO) = QqW\0), where Qq is the detuning from the magnetoexci- 
ton level. On the other hand, as we will see later, the non/mear optical response is dominated 
by the exciton-exciton interactions. 
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Let us now turn to the two-exciton states which govern the FWM signal. Only states 
with zero total momentum contribute to the optical response. The normalized two-exciton 
basis with zero total momentum can be chosen as 

¥^\v,,r2;r[,r',) = ^^ir^,r2)^^^{r[,v',). (11) 

Correspondingly, the states |p, — p) = d^p\0) x c/lp|0) form a complete basis set of the 
zero-momentum two-exciton Hilbert subspace. Using Eq. (|^), one can easily expand any 
two-exciton state in this basis. For example, 

44|o) = 4|o) X 4|o) - 1^410) X rfLpio). (12) 

p 

The matrix elements of the Coulomb potential between four photoexcited particles, two 
electrons and two holes, evaluated in the two-exciton basis (0), have the form 

Ke(p, q) = v\p-^\e-^^-''^'''/'-'^^^'^y^'\ (13) 
VhHip, q) = ^;|p-q|e-(P-'^)'''/2+^(PX'^)■^'^ (14) 
Vehip, q) = -2i;|p_q|e-(P-'^)'''/2 + 26,5^^, (15) 

with Ep given by Eq. (0) . The two-exciton matrix elements of the Hamiltonian H are then 
given by 

H{p, q) = 2fi/pq - V{p, q)/L2. (16) 

The first term is simply the energy of two non-interacting excitons (in the rotating frame). 
The second term represents the exciton-exciton interaction potential 

F(p, q) = 4t;|p_q| sin2[(p x q) ■ z 1^2] e-^^~"^^"'" + Va{v, q) = ^-^(p, q) + Va{v, q). (17) 

The potential Va(p, q) describes the contribution to the exciton-exciton interactions in the 
quasi-2D case that comes from the lifting of the electron-hole symmetry in a strong magnetic 
field by, e.g., the differing band offsets, the lateral confinement, and the disorder.Sa Such 
an asymmetry can be strong in typical quantum wells and plays an important role in the 
nonlinear optical response in QW's.@§ 

Note here that the symmetric part of the exciton-exciton potential, KslPjq)? and the 
exciton dispersion are not independent, but are related as follows: 

2^p-2eo = / (^^s(P.q)- (18) 

This relation reflects the two possible ways of arranging two electrons and two holes into 
two excitons. 

Note also that Vs(p, q) vanishes for zero momenta. This reflects the fact that, due 
to the electron and hole symmetry in the ideal 2D case, the net interaction between two 
magnetoexcitons vanishes at small e-h separation (which corresponds to small momenta). 
In contrast, the potential VA(p,q), originating from the breaking of the e-h symmetry in a 
is finite for small momenta and leads to the main interaction-induced contribution 
to the FWM signal.0ii 
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III. THIRD-ORDER OPTICAL RESPONSE OF 2D MAGNETOEXCITONS 



As discussed in the Introduction, the calculation of the magnetoexciton nonlinear opti- 
cal response in QW's is simplified considerably due to the "freezing" of the kinetic energy 
of electrons and holes by strong perpendicular magnetic field. Using the canonical trans- 
formation method,Elil outlined in Appendix ^ we have obtained a simple expression for 
the FWM polarization Pxx that comes from the exciton-exciton interactions. Refering the 
reader to Appendix O for the derivation, the result reads 



Pxx{t,T) =tfi^£i9{t + T)e- 



(19) 



Here, Si{t) is the electric field of the probe (i = 1) or pump (i = 2), /i is the interband 
dipole matrix element, Uq is the central optical frequency, Qq is the pump detuning from the 
magnetoexciton level, F is the magnetoexciton dephasing rate, and 6{t) is the step-function. 
Since the FWM polarization along the direction 2k2 — ki is linear in the probe optical field, 
we assume, for simplicity, a delta-function probe pulse Si{t) = e^^'^°'^£i6(t + r) and a pump 
pulse of finite duration to centered at t = 0. 

As can be seen from Eq. (|19]), the interaction-induced FWM polarization is directly 
proportional to the two-exciton amplitude x(t). In Appendix ^ we show that, in the case 
of magnetoexcitons in the LLL, x(t) satisfies a simple equation: 

tdtxit) - 2nox{t) = V'^p^t) + F{t), (20) 

where p{t) is the linear pump polarization, determined by Eq. ( P19|) , and V^^ is the Hartree- 
Fock interaction given by Eq. ( |B25| ). The last term in Eq. ( pO] ) describes the effect Of 
four-particle correlations and has the form (see Appendix |BD 

F{t) = - fdciVf^w.it), (21) 

IT J 

where Wp{t) are the coefficients of the expansion of the two-exciton wave-function in the 
basis ([lT|). As shown in Appendix , they satisfy the Wannier-like equation 

tdtWpit) - 2n,w,{t) = - j ^1^(P, ^)w,{t) + V(t), (22) 

where V^^ is the Hartree-Fock potential whose explicit expression is given by Eq. ( P22|) . The 
amplitude x(^) is related to Wp{t) simply as x(^) = WQ{t) — J2p'^p{'t) (see Appendix P). 
The source (last) term in the rhs of Eq. ( P2[ ) describes the interaction of two pump-induced 
polarizations, p^(t), via the Hartree-Fock potential V^^^ , while the first term describes the 
effect of the exciton-exciton correlations. 

The Hartree-Fock approximation corresponds to setting F{t) = in Eq. (^), which 
yields 

;^^^(t) = -iV""^ f rft'e-2^^o(*-*')p2(t'), (23) 

J —oo 

with the pump polarization p{t) given by the solution of Eq. (|B19| ), 
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(24) 



The FWM polarization within the Hartree-Fock approximation is then obtained by substi- 
tuting x^^it) iiito Eq. (0). The relation between the above formalism and that of Refs. 
P],|55|,pB| is discussed in Appendix D. 



IV. ANALYTIC SOLUTION IN STRONG MAGNETIC FIELD 

The role of the exciton-exciton correlations in the FWM polarization is described by 
the function F{t), given by Eq. (|2l|). In order to obtain F{t), it is necessary to find the 
two-exciton amplitude Wp{t) by solving Eq. (|2^). In the general case, this can only be done 
numerically, and the corresponding results are presented in the next section. In this section, 
we derive an analytical expression for Wp{t) in the case of strong magnetic field. 

As discussed above, the interactions between QW magnetoexcitons, confined to the 
LLL, contribute to the nonlinear optical polarization only in the presence of electron-hole 
asymmetry.@'0 The corresponding potential Va(p, q) is short-ranged, and its specific form 
depends on the sample.S In order to proceed further, we assume the following form for the 
s-wave component of this potential: 

^A(p,q) = V:4e-(^'+'^')'^'/^ (25) 

where a is the potential range. The potential strength is characterized by the energy scale 



With such Va(p, q), the Hartree-Fock parameters, determined by Eqs. ( |B22|) and (|B25D , 
have a simple form, 

^r = «o(l-^)e--^«V2. (27) 

and 

Turning to Eq. (p2D, We observe that a considerable simplification occurs for sufficiently 
strong magnetic field such that I < a. In this case, the characteristic momenta of the 
excitons scattered by the potential ( P^D are small: p ~ a^^ < l^^. Since for small momenta 
the magnetoexciton size is also small (see Section II), the symmetric part of the exciton- 
exciton potential, Vsip, q), is suppressed. Indeed, for the characteristic momenta p, g ~ 1/a, 
we get from Eq. ([T7|) that Vs(p, q) ~ e^a(//a)^. Therefore, V5(p, q) can be neglected, as 
compared to VA(p,q), under the condition 




Wo>^o(-) , (29) 
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where Eq ~ c^/kZ is the magnetoexciton binding energy [see Eq. (H)]. Note that the above 
condition can be met even for weak asymmetry, uq <^ Eq, provided that the magnetic field 
is sufficiently strong. Thus, for / < a, we can replace V^(p,q) by VA(p,q) in the first term 
in the rhs of Eq. (^21). Then, by using Eqs. (|2l|), (|25|) , and (^7]), this term can be expressed 
via F(t) as 



12 g-p2a2/2 



(27r)2 ^^^'^^ ' -K^P/a^ - 
Eq. ( p2|) can be solved by Fourier transform: 

-p2a2/2 



- / rfql^f ^^,(t) = -F{t) f——^. (30) 



Fiu) 



+ 



(31) 



where p*^^^ {uj) is the Fourier transform of p'^{t). To find -F(co'), we multiply the above equation 
by V^^ and take the sum over p. This gives 



F{oj) 



P Fjoo) 
'nP/a^- 



P 



-P^o?/2yHF 



dp 



HF\2 



UJ 



2a 



yielding 



where 



Fiio) 



[UJ 



Qiuj) 



l + Qiuj) 



Qiuj) 



Va 



dp- 



UJ 



2a 



Using Eq. (^3]), the amplitude xi.^) can then be easily obtained from Eq. 

yHF p(2)(^) 



as 



UJ — 2a 



S{uj) 



where 



S{UJ) 



l + Q{uj) 



(32) 



(33) 



(34) 



(35) 



(36) 



Equations ( p4j36|) , which are the main result of this section, provide an analytic expression 
for the two-exciton amplitude xi^)-, which determines the FWM polarization (0). The 
result (|35| ) has the form of a product of the Hartree-Fock result, Eq. (|23|) , and the frequency- 
dependent factor S{uj): x{'^) = X^^ {'^)S{uj). In the time domain, x(t) has the simple 
form 



/oo 
dt'S{t - t')x''^{t') 
-oo 



(37) 



with x^^it) given by Eq. (23). As can be seen, S{t) represents a memory kernel whose time 
dependence is governed by the exciton-exciton correlations. In the absence of correlations, 
corresponding to Q{uj) = in Eq. (p6D, the kernel is instantaneous, S(t) = 6{t), and we 
recover the Hartree-Fock result. 
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V. MEMORY FUNCTION 



In this Section we analyze the properties of the memory function ^(ti;), given by Eq. 
(^). For sufficiently strong magnetic fields such that / < a, the main contribution to the 
integral in Eq. (0) comes from small momenta q ~ < . For such momenta, the 
exciton dispersion is quadratic, fig ^ + g^/2M, where M is magnetoexciton effective 
mass [see Eq. (|])]. We then get from Eq. (p4|) 



oo 



~E/D 



Q{uj)=VaU dE (38) 

JO uo — Z\Iq — ibj 



where M = M/27i = -^{EqI'^) ^ is the 2D magnetoexciton density of states and 

2 

2Ma2 ^ 16 I a 



B~ = 2^-l (39) 



plays the role of the cutoff energy. 

For large frequencies, Icj — 2^0! ^ D, the integrand in Eq. (^) can be expanded in terms 
of {lj — 2Qq)~^ yielding 

For small frequencies, \uj — 2Qq\ ^ D, the integral in Eq. ( PBD diverges at the lower limit 
and can be estimated as 



Q{uj) ~ -Gin 

where 



2D 



uo — 2f2o 



-i-KGQiuj -2^^), (41) 



is a dimensionless parameter characterizing the strength of the exciton-exciton interaction. 
As can be seen from Eq. (|4^) , the magnitude of G can be tuned by varying the magnetic field: 
G increases with decreasing the magnetic length /. Note also that the logarithmic frequency 
dependence of Q{lS) is similar to that of the exciton-exciton scattering matrix.i30 

In the following, we will distinguish between attractive and repulsive exciton-exciton 
interactions. 



A. Attractive interaction 

We first consider the case of attractive exciton-exciton interactions, corresponding to 
G > 0. For Lo < 2Qo, the imaginary part of Q{uj) vanishes and the memory function S{lj) 
exhibits a biexciton pole. In the case of weak exciton-exciton interaction (weak asymmetry), 
the corresponding frequency uJh can be determined by using the low-frequency asymptotics 
(13) of Qioo): 1 - Gin = 0. This gives 
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(43) 



where is the biexciton binding energy. Note that the above expression holds for G <^ 1 
corresponding to small biexciton energy <^ D. The biexciton contribution to S{uj) can 



be obtained by expanding Eq. (|41| ) in the vicinity of uj^: 



or, in the time domain, 



Si(t) = !^9(i)e-l«>-2««)'-i.', (45) 



where 9{t) is the step function and we have included the biexciton width 7f, that mainly 
comes from electron-phonon processes not incorporated in the Hamiltonian H considered 
here. 

In the time domain, the biexciton resonance leads to oscillations of the memory kernel 
with a period determined by the biexciton binding energy (see Eq. (^)). The memory effects 
decay on a time scale determined by the biexciton dephasing time, which, in a magnetic 
field, can be considerably longer than the exciton dephasing timc^'il Near the biexciton 
resonance, the amplitude xi.^) [see Eq. (|35D] that determines the FWM signal, takes the 
form 

yHF {2)t \ 

M^) = \r,\l\- . (46) 

uj - 2\Iq + Eb + ijb 

and exhibits a biexciton pole as a function of frequency. 

The memory function S{uj) also exhibits a second resonance centered at the frequency 
ujsc = 2f2o + Efy, which describes the continuum band of exciton-exciton scattering states. 
The width of this broad continuum resonance is determined by the exciton-exciton interac- 
tions. Indeed, for u > 2Qq, the function Q{uj), that describes the exciton-exciton scattering, 
develops an imaginary part proportional to the interaction strength [see Eq. (|41|)]. For 
Eb > 7fe, the scattering band and biexciton resonance are well separated from each other 
and lead to distinct dynamics, while for Eb < 76 they merge into a single asymmetric peak. 

The above results hold in the case of weak exciton-exciton scattering, G -C 1. With 
increasing magnetic field, the parameter G, Eq. (^2p, that characterizes the strength of 
the exciton-exciton interaction increases, and so does the biexciton binding energy Eb- For 
sufficiently strong magnetic fields such that Eb ^ D, corresponding to G 3> 1, one can use 
the high frequency asymptotic expansion of Q{uj), Eq. (^0]), which when substituted into 



Eqs. (|3g) and (|35D gives 



Note that for G ~ 1, we have Eb ^ Uq, and Eqs. (|47| ) and ( ^q ) match. 

It is important to note that Eq. ( ^Tj) becomes exact in the strong field limit, / -C a. 
Indeed, in this case, the typical magnetoexciton momenta are small, q ~ -C so that 
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one can disregard the magnetoexciton dispersion in Eq. (|3^ and replace Qq by Qq, which 
leads again to Eq. (^OD for Q^u) and, therefore, to Eq. (^7|) for x{^)- Iii other words, for 
G S> 1, the exciton-exciton scattering is suppressed and the memory effects are dominated 
by the biexciton resonance. The simple expression ( ^7|) can be then viewed as a "biexciton 
pole approximation" for the two-exciton amplitude x{^)- fact, Eq. (^) is equivalent to 
the average polarization modeflffiS'i (APM), as we show in Appendix 0. 

In Fig. (0) we plot the memory function, obtained by evaluating Eqs. (^) and ( |5^ 
with the full exciton dispersion (|^). As we can see, for attractive iteraction, ^(ti;) is mainly 
dominated by the biexciton pole. 



B. Repulsive interaction 

Let us now turn to the case of repulsive exciton-exciton interactions, G < 0. In this case, 
the memory function describes a continuum band of scattering states above the two-exciton 
energy. We start with the case of strong interactions \G\ ^ 1. Here we can use again the 
high frequency asymptotics ( PP] ) for ReQ^u) (valid for u — 2Qo ^ D). We see from Eq. 
that 5'(ci;) exhibits a resonance at u = 2Qq + Imo]. In contrast to the attractive case however, 
in this frequency range Qi^uj) develops imaginary part, which can be deduced form Eq. (|38D 
as lmQ{u) = —i7cGe~^'^~''^^'^^^'^^9{uj — 2i7o)- Near the resonance, we obtain from Eqs. ( |36[ ) 
and Eq. (^) that 

XM = 1, ^1 71. , 7sc = ™oGe-l^l (48) 

uj — 2\Iq — \uq\ + 



I sc 



The above expression describes a resonance made up from the continuum of scattering states. 
Eq. (^Hj) is valid for \uq\/2D = \G\ ^ 1, indicating that in this case the resonance is sharp. 
Comparing to Eq. (^Sp, we can view this continuum resonance as an antibound state. The 
crucial difference, however, is that here the resonance width is determined by the exciton- 
exciton interactions rather than homogeneous broadening due to phonons. The narrow width 
of the resonance comes from the fact that for |G| ^ 1, the momentum exchange processes 
between excitons are suppressed. Note here that the APM ( ^Tf ) is equivalent to Eq. (^Sf) but 
with jsc introduced phenomenologically. 

With decreasing the lineshape of the memory function changes. This is illustrated in 
Fig. (^. The sharp Lorentzian peak transforms into a broad asymmetric band centered at 
lower frequencies. Importantly, S{uj) develops a cusp at a; = 2f2o- This cusp is manifestation 
of the logarithmic singularity in the low-frequency asymptotics of Q{uj) [see Eq. (^)]. In 
the next section we will see that this low-frequency behavior of the memory function has a 
dramatic effect on the FWM polarization for long time delays. 



VI. NUMERICAL RESULTS AND DISCUSSION 

In this section, we present the results of numerical calculations of the third-order FWM 
polarization and compare them to the analytical solution of Sec. |V| and |V[ The calculations 
below also incorporate the Pauli blocking contribution, PpB{t,T), to the total FWM polar- 
ization P(t, r) = Pxxit, t) + Ppsit, 't), which is given in Appendix 0. We use the standard 
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conventions according to which the positive time-delay, r > 0, corresponds to probe pulse 
arriving before the pump pulse. We will be interested mainly in the time-integrated and 
spectrally-resolved FWM signals, defined as 



TI-FWM 



/oo _ 
dt\P{t,T)\' 
-oo 



(49) 



and 



SR-FWM = |P(u;,r)p, 



(50) 



respectively, where P{u!,t) is the Fourier transform of P{t,T). 



A. Attractive interaction 



We start with the case of attractive exciton-exciton interactions, Uq > 0. In Fig. § we 
show the results of our numerical calculation of the TI-FWM signal for a strong magnetic 
field, a/l = 6.0, and for both strong and weak e-h asymmetry. In order to be consistent with 
the recent experimental data|iS we have chosen the biexciton homo geneous broadening, 
76, to be smaller than that of the exciton, F. In both cases, the calculations yield a strong 
TI-FWM signal at negative time delays. For strong asymmetry (large uq), the biexciton 
oscillator strength in the memory function S{u!) is enhanced (see Fig. ([l|)), and the TI-FWM 
signal exhibits biexciton oscillations. This behavior can be easily deduced from the analytic 
results of the previous section. Indeed, for 5-function pump pulse centered at t = 0, the 
pump polarization ( ^If ) is given hj p{t) = —ifiS20{t)e~^^°'^~^'^, and the biexciton contribution 
to the amplitude x(^) [Eq. (P6|)] takes the simple form 



The first term in the rhs leads to the oscillations in the TI-FWM signal with the period Tf, = 
271 /Eb, and to the overall exponential decay of the TI-FWM signal with the characteristic 
time at long negative time delays. Such a behavior is also consistent with the calculations 
in Ref. 

In Fig. |, we plot the time evolution of the corresponding SR-FWM signal, which shows 
both an exciton and a (weaker) biexciton peak. A distinguishing feature of the spectra 
is that the height of the exciton peak oscillates as a function of time delay with a period 
determined by the biexciton binding energy, Tb = 2ti/ Eb. This behavior is easily reproduced 
by using Eq. (^) (see below). 

As mentioned in the previous section, in the strong magnetic field case l/a -C 1, the 
exciton-exciton scattering is suppressed and the FWM polarization is governed by the biex- 
citonic effects. In Fig. ^, we show the TI-FWM signal for weaker magnetic field such that 
a/l = 3.0 and for two different pump durations, both of which are shorter than the exciton 
dephasing time, to < T"^- In Fig- @(a), the pulse duration was chosen to be shorter than 
the inverse biexciton energy, < E^^. In this case, the TI-FWM curves are qualitatively 
similar to those in Fig. 0(a); in particular, they exhibit an overall exponential decay at long 
negative time delays together with biexciton oscillations as Uq increases. In contrast, for 




(51) 
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the longer pulse duration to > -E'jT^; the change in the shape of the TI-FWM signal with 
increasing Uq is different [see Fig. ^(b)]. 

To understand this behavior, we note that, for weak asymmetry uq/Eq <C 1, we have 
to < -E^^, and the decay of the TI-FWM signal is similar to that for the shorter pulse 
duration [compare Figs, ^(b) and |(a)]. On the other hand, for larger uq we have to > ^ 
and the pump pulse tuned at the magnetoexciton level, fio = 0, does not directly excite 
the biexciton. In this case the decay at intermediate times becomes non-exponential , while 
the amplitude of the oscillations is significantly enhanced [see Fig. |^(b)]. The reason for 
such a behavior is that, for short time delays, the signal decays exponentially with the 
Hartree-Fock decay timeO (4r)~^ [see Fig. |5(b)], while for longer time delays, the build-up 
of the biexciton correlations leads to an exponential decay with a different characteristic 
time, (27f,)~^, determined by the biexciton width. The crossover between these two regimes 
leads to an apparent nonexponential decay at intermediate time delays [see Fig. |3(b)]. 

The sensitivity to the detuning of the time-evolution of the FWM polarization becomes 
even more evident in the SR-FWM signal. In Figs. § and ^ we plot the corresponding 
spectra versus time delay for pump tuned to the exciton and biexciton levels, respectively. 
For VLq = 0, the exciton peak decays exponentially for positive time delays, r > 0, with 
characteristic time (2r)^^, while the biexciton peak is suppressed (see Fig. |^). For pump 
tuned at the biexciton, 2Qo = Ef,, the situation changes: for r > the decay of the 
exciton peak follows the pump pulse, while the biexciton peak decays exponentially with a 
characteristic time (2r)~^ (see Fig. |^). For negative time delays, the exciton peak oscillates 
with the period 2n/Eh (compare to Fig. ^). 

To understand this behavior, consider the biexciton contribution to the frequency- 
dependent polarization Pxx{^^,t). The latter can be presented as a sum of two terms 
P^^\uj,t) + P/^''(c<j, r), corresponding to the first and second terms in Eq. ([19|) . The first 
term can be presented as 

/^"(.■r)^- '''^""!?!;^''-' . (52) 

UJ + ZL 

where uj = uj — Eg + Eq is the frequency measured from the magnetoexciton level (we used 
here that Qq = Eg — Eq — uq). Equation (|5^ ) describes the peak at the exciton energy, 
uj = 0. The peak height depends on the time delay via Xb{—T)- The latter, for short pulse 
duration, can be approximated by Eq. ([51|), leading, for r < 0, to the oscillations of the 
exciton peak amplitude with the period 2TT/Eb determined by the biexciton binding energy 
(see Figs. ^ and |). Note that, since Xb(t) = for times t < prior to the onset of the pump 
pulse [see Eq. (|5lD ], the above contribution vanishes for r > 0. Therefore, for positive time 
delays, the SR-FWM signal at the exciton frequency is primarily determined by the Pauli 
blocking term. 

The second term in Eq. (Il9[) has the form 



r) = ,^£,e-<-«.>' r '^^^pi^. (53) 

j-oo 2tt uj + 2\lo — uj' + tl 



Using Eq. (RBI), the integrand of Eq. (|55D can be decomposed as 



UJ + 2^0 - u' + iV UJ + Eb + iiV + -fb) 



1 1 



uj' - 2no + Eb + i-fb UJ + 21^0 -uj' + iT 



(54) 
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yielding 



~(2) _ /i2^iG-ie~*('^+^o)^ duj' 



+ 



uj' - 21^0 + Eb + i-fb + 2^0 - a;' + 



(55) 



The prefactor in the rhs of Eq. has a resonance at the biexciton frequency, uj = —Eb, 
with width F + 7f,. The amphtude of the biexcitonic peak is determined by the integral in 
Eq. For negative delays, the main contribution to the above integral comes from the 

first term in the integrand, whose denominator has a pole at u' = 2Qq — E^. For resonant 
pump excitation, VLq = 0, and pulse duration to < T \ the frequency width of p*^^^ (cj') in 
the integrand of Eq. (^) is of the order of 2F. In this case, the contribution of the pole is 
suppressed for Eb > F, and therefore the biexciton peak is diminished (see Fig. For off- 
resonant pump tuned at the biexciton, 2^0 — Eb, and for pulse duration to > 2/Eb ^ i^o\ 
the time dependence of p{t) is mainly determined by that of the pulse itself, so that the 
frequency width of p^'^^tu') is ~ t^^. In this case the pole dominates the integral, leading to 
a strong biexciton peak in the SR-FWM signal, which decays on a time scale (27b) ~^ (see 
Fig. 0). 

For positive delays, r > 0, the main contribution to the biexciton term P^; (cu, r) comes 



from the second term of the integrand in Eq. (^5]), which exhibits a pole at u' = Cu + 2^0. 
Since O ~ —Eb due to the prefactor in Eq. (|53|), for resonant excitation, = 0, the biexciton 
peak is diminished for > F (see Fig. For off- resonant excitation, 2f2o = Eb, this pole 
dominates the integral leading to a strong SR-FWM signal at the biexciton energy (see Fig. 
1^, which decays on a time scale (2r)~^ [see Eq. (pSj)]. Note that, for pump excitation at the 
biexciton frequency, 2f2o = Eb, the SR-FWM signal at the exciton frequency is strong for 
negative time delays due to the exciton-exciton correlations [see Eq. (p^j, but diminishes 
after the pump pulse is gone due to the suppression of both the Pauli blocking contribution 
and P^^^ [see Eq. (ET)]. 



It should be noted that for all the above parameter values, the analytical curves calculated 
from Eqs. (|5^P^), which were derived in the previous section for I < a, are practically 
indistiguishable from the exact numerical calculations in Figs. ^-0. In fact, Eqs. (|3^ - |36D 
provide a very good approximation even for a ~ /. This is demonstrated in Fig. H, where we 
present the results of numerical calculations of the TI-FWM signal at a// = 1.5 together with 



the analytic solution ( ^^ - |36|) for different values of the biexciton width 7^. For comparison, 
we also plot the TI-FWM signal derived using the APM {^7\). Note that, for the parameters 
used in Fig. |^, Uq/D ~ 5.0 and we are therefore in a regime where Eq. (^) is meaningful. We 
see that our analytic solution (P^D, which includes the effects of the exciton-exciton scattering 
processes, neglected by the APM, is in a much better agreement with the numerical results 
even for strong e-h asymmetry. In particular, the APM (^) fails to provide the correct 
period and amplitude of the biexciton oscillations for both 7^ -C F or 7^ = 2r. In addition, 
the APM gives a considerably weaker magnitude of the FWM signal. 



B. Repulsive interaction 

Let us now turn to repulsive exciton-exciton interactions, Uq < 0. In Fig. ^ we show the 
TI-FWM signal in a strong magnetic field, a/l = 3.0, for several values of the parameter 
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G [Eq. (|4^)] characterizing the strength of the exciton-exciton interactions. For large \G\, 
the signal decays exponentially with a characteristic time that increases with In this 
regime, the memory function S{uj) is dominated by an antibound state well separated from 
the exciton resonance (see Fig. |^). A narrow antibound resonance leads to oscillations in the 
TI-FWM signal, similar to biexciton oscillations. Note that antibound states were recently 
observed experimentally in the FWM signal of QW's at zero magnetic field.0 The decay 
time is simply the inverse antiresonance width, given by Eq. (^). Importantly, the latter is 
determined by the exciton-exciton interactions, consistent with Fig. ^. 

As \G\ decreases, we observe a transition from an exponential regime for initial time 
delays, where the antiresonance dominates the decay, to a nonexponential regime for long 
negative time delays. For small the latter regime extends down to times ~ 2/r; for 
longer times, the TI-FWM signal curves up with increasing |r| (see Fig. |^). 

The long nonexponential tail of the TI-FWM signal has its origin in the low-frequency 
behavior of the memory function S{lj). As can be seen from Eq. (|r^), the shape of the signal 
for long negative delays is determined by the long-time asymptotics of the amplitude xi^)- 
The latter can be derived by using the analytic expression ( pSf ) for x{^)- For sufficiently large 
t, we can substitute into Eq. (|35D the low-frequency asymptotics of the memory function 
S{uj) [see Eq. (|T|)], 

Si^) ^ —7:7. on ^ . 7777- (56) 



1 + +^^\G\e{uJ-2^lo) 
As shown in Appendix 0, for short pump duration this gives 

^(*) - ^iW^' '''' 

with D given by Eq. (^). The above result is valid for long times t ^ D^^. The slow 
logarithmic decay of the amplitude xit) reflects the singular behavior of >S'(ti;) at small 



frequencies. Note that both the real and the imaginary part of the memory function (|56D 
contribute to the asymptotics ( |57D (see Appendix]^). 

Using Eqs. (^) and (|19]), the TI-FWM signal (|4^) for negative time delay takes the form 

TI-FWM = ^i^v^^^ / t/te-2r(*-l-l) ^ _^ . (58) 

412 _/|^| + |G|ln2D|r| l + \G\\n2Dt) ^ ' 

The main contribution to the integral comes from times (t — |r|) ^ F^^. For |r| ^ F"^, the 
second term in the rhs of Eq. (|58|) can be expanded in {t — \t\)/\t\, and we finally obtain 

TI-FWM oc —— (59) 

The above expression describes the long tail of the TI-FWM signal for negative time delays 
in Fig. 1^. Note that Eq. (^) is valid for |r| 3> F"^, so that with increasing F, the crossover 
to the asymptotic behavior ( p9D occurs at earlier time delays, as can be seen in Fig. |^(b). 

The behavior of the TI-FWM signal on various time scales can be tuned by changing the 
magnetic field strength. For strong magnetic fields / -C a, the momentum exchange processes 
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between excitons are suppressed and the TI-FWM signal is dominated by the antibound state 
resonance leading to the exponential decay. With decreasing magnetic field, the energy scale 
D which determines the characteristic exciton momenta q ~ (2MZ))^/^, increases [see Eq. 
(^)], and the crossover to the above nonexponential regime occurs at earlier times |r| > 
[see Figs. |10|(a) and|TT](a)]. At the same time, the shape of the TI-FWM signal for |r| > D^^ 
depends strongly on G. For |G| ~ 1, the decay is described by the equation 

TI-FWM oc 'p J m 

For IGI ^ 1, the situation is more intricate. In this case, the asymptotics ( pD| ) applies for 
very long time delays, |r| > D^^e^ ^ D^^, while for D^^ < |r| < D^^e^ the decay is 
quadratic. 

The behavior of the TI-FWM signal at initial times, |r| ^ D^^, is mainly determined 
by the imaginary part of S{uj). At small the memory function is a smooth function 
of frequency for > D (see Fig. ||), and, correspondingly, S(t) is nearly instantaneous, so 
that the TI-FWM signal shows a Hartree-Fock-like decay (see curves with \G\ = 0.24 in Fig. 

At IGI ~ 1, ImS'(co') represents a broad asymmetric peak centered at E^c ~ Uq with the 
width '-fsc ~ -D, so that TI-FWM signal exhibits weak oscillations and overall exponential 
decay on time scale |r| ~ D~^. For long time delays, |r| ^ D~^, the low-frequency singular 
behavior of the memory functions becomes important and the crossover to the asymptotic 
regime Eq. (^) occurs (see curves with \G\ = 1.2 in Fig. |^). For larger the crossover 
occurs at even longer times, |r| ^ D^^ (see curves with \G\ = 3.6 in Fig. 

In Figs. |10| and |ll] we compare the numerically calculated TI-FWM signal with the 
analytic solution (|3^ - |36D and with the APM (^). In the latter case, the width of the band 
of exciton-exciton scattering states, 7sc, must be introduced phenomenologically by fitting 
the memory function S{uj) with a Lorentzian [see Figs. |lO|(b) and pA|(b)]. Even though 
such an approximation describes well the overall lineshape of the memory function, the 
latter deviates from a Lorentzian in the low frequency domain. These discrepancies have 
a profound effect on the TI-FWM signal. Namely, even though in the initial exponential 
regime the APM and numerical curves are in a reasonable agreement with each other, the 
APM completely fails to reproduce the long nonexponential tail of the TI-FWM signal. The 
latter originates from the singular low-frequency behavior of S{uj) [see Eq. (pB])]. On the 
other hand, the analytic solution (|3^ - |36D describes accurately the long nonexponential tail 
of the TI-FWM signal as well as the crossover between the two regimes [see Figs. p!0|(a) and 
0(a)]. 



VII. CONCLUSIONS 

In summary, we studied the role of the exciton-exciton interactions in the ultrafast non- 
linear optical spectroscopy of semiconductor quantum wells in a perpendicular magnetic 
field. In the case of attractive exciton-exciton interactions, we found that the biexcitonic 
effects dominate. For repulsive interaction, we have shown that the time-integrated four- 
wave-mixing signal exhibits a long nonexponential tail for negative time delays. We traced 
the origin of this tail to the low-frequency of the memory kernel for the two-exciton wave- 
function, for which we derive an explicit analytical expression. 
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Our analytical solution, given by Eqs. (|3^ - |36D , was derived here in the case of strong 
magnetic field, where the kinetic energy of the electrons and holes is frozen and the exciton- 
exciton interactions play a dominant role in the third-order optical response. We believe, 
however, that Eqs. (|3^ - |36D can serve as an analytical model also for the zero magnetic 
field case, for example, when the dynamics is dominated by interactions of Is excitons. In 
particular, our solution reproduces, as a limiting case, the well-known average polarization 
model,00ii'0ii which has been used in the case of zero magnetic field. Unlike the APM, 
however, the model (|34| - p6D takes into account nonperturbatively the correlation effects due 
to the continuum of exciton-exciton scattering states. 
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APPENDIX A: 



In this appendix, we present a general expression for the third-order FWM polarization 
derived using the canonical-transformation formalism.ElEBEl The total Hamiltonian of the 
system is Hfot = H + Hp + Hs, where H is the two-band semiconductor Hamiltonian (consid- 
ered in the rotating frame from now on) and describe the coupling to the optical fields 
(a = 1, 2 for probe and pump, respectively). 



Ha = -M«(t)[/^ + H.C. 



(Al) 



with 



M„(t) 



(A2) 



Here f/^ is the optical transition operator, aj^ and b^_^ are the creation operators for electron 
and hole respectively with momentum k; Sa{t) are the amplitudes of the optical fields 
propagating in the probe (ki) and pump (k2) directions with central frequency ujq, and fi is 
the dipole matrix element. 

The nonlinear optical polarization is given by P(t) = fi{'^(t)\U\'^(t)) , where the state 
|\E'(t)) satisfies the Schrodinger equation with Hamiltonian Htot- The third-order polarization 
propagating along the FWM direction 2k2— ki can be obtained by expanding the state {"^{t)) 
to the first order in the probe and second order in the^ump field. The FWM polarization 
then takes the form PpwMit) = e*(2k2-ki)-rp^^) withH 



Pit) 



iuJot 



dt'£*{t') 



FWM 



(Olo)-(t-O 



(A3) 



The above expression has a form similar to the linear polarization due to the probe optical 
field, with the important difference that here the single e-h pair state UpyyM{'t)\^) depends 
on the pump optical field. 
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uIwaM = UW^\0) - VW\0), (A4) 

where |0) is the ground state of H. Here the states V^(t)\0) and Vr^(t)|0) are the first and 
second order terms of the expansion of the state that evolves from the semiconductor 

ground state under photoexcitation by the pump alone, 

oc |0) - e'''^"pt(t)|0) + e^'^' ''W\t)\0). (A5) 

The single e-h pair state 'P'''(t)|0) and the two e-h pair state iy''(i)|0) excited by the pump 
optical field satisfy the equations 

idtV\m = HV\tm+^£2{t)U^\Q), (A6) 

and 

idtW\t)\Q) = HW\t)\Q) + ii82{t)U^V\t)\Q), (A7) 

respectively.S'El 

The decomposition ( |A^ ) has a straightforward interpretation. The first term in the rhs 
describes the FWM signal coming from the excitation of two e-h pairs by the pump followed 
by the de-excitation of one pair by the probe. The second term describes the excitation of an 
e-h pair by the pump, followed by de-excitation by the probe and then excitation of an e-h 
pair by the pump again. The first term in Eq. ([A4|) describes the two-exciton contribution 
to the FWM signal, while the second term describes the one-exciton contribution. Note 
that, in a two- level system, only the second term in Eq. ([A4|) is present. 



APPENDIX B: 

Starting from the general formulae in Appendix ^ we derive here a simple expression 
for the FWM signal coming from the exciton-exciton interactions. The Pauli blocking con- 
tribution is outlined in Appendix |^. 

Equation (|A7|) describes the time evolution of the two-exciton state from an initial state 
of two non-interacting excitons photoexcited by the pump. It is useful to separate out 
the contribution Wd(t)|0) due to non-interacting excitons with zero momentum by writing 
W\t)\Q) = W^{t)\0) + W4x(^)|0), where the second term comes from the exciton-exciton 
interactions. In the following we assume that the photoexcited electrons and holes are 
confined to their respective LLL's, so that the energy of two noninteracting magnetoexcitons 
with zero momentum is simply 2Qo, where Qo is the pump detuning from the zero- momentum 
exciton energy (we work in the rotating frame). 

The time evolution of the non-interacting exciton state is determined by the equation 

tdtW^{t)\0) = 2QoW^{t)\0) + fiS2{t)U^V\t)\0). (Bl) 

After subtracting Eq. ( pi| ) from Eq. ([A7| ), we obtain that the exciton-exciton interactions 
are described by the time-dependent Schrodinger-like equation 

zdtWj,^{t)\0) = HWj,^{t)\0) + {H- 2no)WMo). (B2) 
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The source term on the rhs is the Hartree-Fock interaction between two zero-momentum 
excitons. The exciton-exciton correlations come from the Coulomb potential in the Hamil- 
tonian H (first term in the rhs of Eq. (P^)). 

Thus, we can present Eq. (^) as a sum of two contributions: 



uIwm\^) = uU^) + u],^\q), 



(B3) 



where 



(B4) 



describes the Pauli blocking effects for noninteracting excitons, and 

Uix{m) = UW],^\Q) (B5) 

comes from the exciton-exciton interactions. The polarization ( [A3[) can then be presented 
as a sum of Pauli blocking and exciton-exciton interaction parts, P{t) = PpBif) + Pxxit), 
where 



2 —iojot 



PpB{t) = ijj e 



dt'£*it' 

o 

dt'Slit') 



(t)PB{t,t') - (l)pB{t\t) 

<l)xx{t,t')-(l)xx{t\t) 



(B6) 



with 



<PPBit,t') = {0\Ue-''''^'-''^uUt')\0), 
0xx(t,t') = (O|t/e-^(*-*')t/.L(t')|O). 



(B7) 



Equations ( P6| ) provide an exact formal expression for the third-order FWM signal. 

Here we focus on the the interact ion- induced FWM polarization, and defer the Pauli 
blocking contribution to Appendix y. Pxxit) is determined by the correlation function 
4>xxit,t'), which satisfies the equation 



tdt<l>xx{t,t') = (0|t///e-^^(*-*')t/j,^(t')|0) 
with the initial condition (pxxitjt) = x{t)j where 

xit) = {o\uuj,Atm = {o\uuwj,Atm. 

Since the state W\0) = N^/^dl\0) is an eigenstate of H [see Eq. 
51 that 

0xx(t,t') = e-^^«(*-*')-^l*-*'lx(t'), 



(B8) 
(B9) 

we obtain from Eq. 

(BIO) 



where F is the exciton homogeneous broadening. Since Eq. ( P10|) is linear in the probe optical 
field, we consider for simplicity a 5-function probe Si{t) = e^*'^"^£^i(5(t + r) Substituting this 
into Eq. (pT]) we obtain 



Pxx{t,T)=zfi^S^9{t + T)e 



-r{t+T)-iuJoit-T) 



-ino{t+T) 



xi- 



jno{t+T) 



xit) 



(Bll) 
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Using Eq. (0), xit) can be expressed in terms of the coefficients of the expansion of the 
state Wxxit)\0) in the two-exciton basis (|TT|) , 



wp{t) = {p,-p\wiAm, 



as follows: 



X{t) = ^0 - ^ E 



Wo - Wq. 



(B12) 



(B13) 



The two terms in the rhs of Eq. (|B13| ) reflect the two possible ways of arranging four particles 
into two excitons. The equation for the amplitude Wp{t) follows from Eq. 



idtwpit) = (p, -p\Hwj,Am + (p, -p\iH - 2no)wMo). 



(B14) 



After projecting out the rhs of the above equation in the two-exciton basis Eq. ( ]TT| ) and 
using Eq. (^) for the matrix elements of H, Eq. ( pi4|) takes the form 



(2 



(27 



(B15) 



where Wp{t) = {p, —p\WQ{t)\0) is the noninteracting two-exciton amplitude. The equation 
for the latter can be obtained straightforwardly from Eq. 



idtw'Jt) - 2(fio - ^r)w''(t) = fiS2{t){p, -p|f/tpt(t)|o). 



(B16) 



In order to evaluate the last term in Eq. ( |B16| ), we use Eq. ( |A6|) for the one-exciton state 
V^{t)\0) to obtain 



(tdt -no + tT){p, -p|f/tpt(i)|o) = ^S2{t){p, -p|f/tf/t|0). 
Using Eqs. (0) and we find that (p, -p\WW\0) = N6po - 1, so that 

'27T 



{p,-p\U^V^{t)\0) = {N6po-l)p{t) 



2 ^(P) - 1 



pit). 



where p{t) is the linear pump-induced polarization, satisfying 



(B17) 



(B18) 



(B19) 



Substituting Eq. ( |B18| ) into Eq. ( |B16| ), and comparing the latter to Eq. ( |B19| ), we obtain 
that 



P\t) 



(B20) 



Substituting the above expression for w^{t) into Eq. ( [BTSD and using the decomposition 
^(P) q) = ^s{p,fi) + Va(p, q) together with the relation ([T8|) , we arrive at the following 
equation for the two-exciton amphtude: 
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zdtWpit) - 2QpWp{t) = -J j^V{p, q)w,{t) + VfVit), (B21) 

where 

According to Eq. ( |B13| ), the amphtude xif) is determined by the p = solution of Eq. ( P21|) 
and by WQ{t) = N~'^ Y^pWp{t). The equation for Wo{t) can be straightforwardly obtained by 
summing Eq. ( P21|) over all momenta p and using the relation (pISl). This gives 

zd.Mt) - 2noMt) = E / (^^a(p, cOw,{t) + i E V (t). (B23) 

Setting p = in Eq. ( [B21D , and subtracting from the latter Eq. ( |B23| ), we obtain a simple 
equation for the amplitude (Eq. (pi3|) ) that determines the FWM polarization, 

zdait) - 2QoXit) = V^V W + Fit), (B24) 

where 

V^^ = Vr-^ldciVf^ (B25) 
is the Hartree-Fock interaction, and 

F(t) = ^ I dqV^''^ w,{t) (B26) 
describes the effects of the exciton-exciton correlations. 



APPENDIX C: 

In this Appendix we derive the Pauli blocking contribution to the FWM signal in strong 
magnetic field, determined by Eq. ( p^ . Using Eq. (|A6|) , the second term in Eq. (|B4|) satisfies 

idtV\t)UP^it)\0) = 2{no-iT)V\t)UP^it)\0)+ fi£2it)[U^UV^it)+V^it)UU^\0) (CI) 

where we used the fact that 'P^(t)|0) is the zero- momentum exciton state with energy Qq. 
Using the relation W = N^/'^d\ [see Eq. (^] and noting that d\dQ = 1 when acting on the 
one-exciton state V^{t)\0), while UU^O) = Ndodl\0) = N\0), we obtain that 

idtV\t)UV\t)\0) = QoV\t)UV\t)\0) + 2Nfi£2{t)V\t)\0). (C2) 

The non- interacting contribution to the two-exciton state, M/d|0), is determined by Eq. ( |B ![ ) . 
After projecting with U on the Ihs we obtain that 

tdtUW^it)\0) = 2noUW^it)\0) + fi£2it)UU^V\t)\0). (C3) 
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The last term can be simplified by using UWV^it)\0) = Ndodlv^ {t)\0) = 2{N - 1)V^ {t)\0) , 
so that 



idtUW^it)\0) = noUW^{t)\0) + 2{N - l)/i^p(t)pt(t)|o). 



(C4) 



Comparing the above equation to Q, we see that UW^{t)\0) = (1 - l/N)V^{t)UV^t)\0), 
and we thus obtain 

uUt)\0) = -^VHt)UVHtm. (C5) 
Using the formal solution of Eq. ( |A6| ) , 

V\t) = -tfi f dt'£2{t') e-^^^(*-*')t/te^^^(t~t')^ (C6) 

and that V'^{t)\Q) = p{t)W\0), where the pump polarization p{t) is given by Eq. (0), we 
obtain 

(j)PB{t,t') = '-^p{t') f rft%(t")(0|t/e-*^(*-*') e-^^(*'-*'Vtei^^(t'-*");7t/t|o). (C7) 

i V J—oo 

Finally, using UW\0) = N\0) and that U^O) is proportional to the magnetoexciton eigen- 
state of H, with eigenvalue Qqj "we obtain that 



(C8) 



The Pauli blocking contribution to the FWM polarization is then calculated by substituting 
the above expression into Eq. (|B6|) . 



APPENDIX D: 

In this Appendix we show the connection between our expression for the FWM signal, 
Eqs. (PUD, (P7|), and ( |B1CI| ), and that of Refs. Q,^,^. The equation for the four- wave- 
mixing polarization Pxxit) = e^'^°'^Pxx{t) can be obtained by taking the time derivative 
with respect to t of the rhs of Eq. (|B6[) , 



(Dl) 



tdtPxxit) = z/i^ /* dt'£*{t') t') - tdt^xxit', t) 

J—oo L 

From Eq. ( piO|) we obtain after some algebra that 

idt(i)xxit,t') = (^]o -«r)0xx(t,t') 

and 



tdt(t)xx{t',t) = (fio - tT)(t)xx{t',t) + 



(D2) 



(D3) 



Substituting Eqs. ( P2|) and ( |D3|) into Eq. (pl| ) and using Eq. (P24| ) we obtain a familiar 
equation for the FWM polarization: 
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idtPxxit) = {Qo - ir)Pxx{t) - fiplit) 



(D4) 



where Piit) is the probe- induced polarization 



(D5) 



The simple (momentum-independent) form of the above equation is due to the dispersionless 
energy spectrum of the electron and hole in the magnetic field. At zero fieldi_a similar equa- 
tion is obtained after averaging out over the spatial degrees of freedom. E^flS The function 

given by Eq. ( P26|) , is determined by the two-exciton amplitude Wp(t), satisfyins 
( [B14D , and describes memory effects due to four-particle correlations. The APMli3£2lE 
is recovered if one uses the simple model (|47| ) for the memory function S{u). 



APPENDIX E: 

In this Appendix, we derive the long-time asymptotics of xi^)- Fo^^ delta-function pump 
centered at t = 0, the pump polarization is simply p(t) = —ifiS20(t)e~^^°^~^^, so that the 
Fourier transform of p'^ (t) is 

uj — 2Uo + 2ii 

Note that for short pump duration, we can write xif) = where x{t) is the 

amplitude for zero detuning; therefore, it is sufficient to evaluate the Fourier transform of 
Eq. (|35|) at Qq = 0. For t ^ F"^, the relevant frequencies are small, -C F, so that 
p(2)(^^~) ~ — (iu£^2)^/2F, and x{t) takes the form 



m = -'-^y^' / — e-*5(^). (E2) 
47ri J-oo uj 

In order to get for long times, t ^ -D^\ we substitute the low-frequency asymptotic 
expression (|56D for the memory function (with f2o = 0), 

m = -^^y''^\m^mi (e3) 

where we separated out the contributions from the real and imaginary parts of S(lo): 



foo duj e"*^*(l + IGllnl^l 
Mt)= -7 \ (E4) 



°° ^ '1 + \G\\n\^\y + n^\G\W{iu) 



J2{t) = -I / — ^ . E5 

Jo u ^ ^ \n\\r. \^-^ ■ 



l + |G|ln|^| +7r2|G|2^(cj) 
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Consider first Ji(t). For t ^ D~^, the characteristic frequencies are small, so that the second 
term in the denominator of the integrand can be neglected as compared to the logarithmic 
term. In this case, the main contribution comes from the imaginary part of the integrand. 



dio sin ujt —in 
~1 + IGllnl^l ~ l + \G\\n2Dt 



Turning to J2{t), we notice that, for t ^ D^^, the main contribution comes from the real 
part of the integrand. After omitting the second term in the denominator and integrating 
by parts, we obtain 



duj sin ut —in 
1 + |G|ln|^| ~ 1 + \G\\n2Dt' 



T / \ / — U/U/OlliU^l^ 1,11 



One can show that the sub leading terms are ~ IGj [^1 + jGI In 2Dtj . Note that the imaginary 
and real parts of S{uj) contribute equally to the long-time asymptotics. Finally, substituting 
Eqs. ([E6| ) and (pT]) into Eq. (p3D, we arrive at Eq. 
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FIGURES 



FIG. 1. The memory function S{uj) for attractive exciton-exciton interaction calculated from 
Eq. (H) with a/l = 3.0 and uq/Eq = 0.05 (solid line), uq/Eq = 0.4 (dotted line), and uq/Eq = 0.8 
(dashed line). The oscillator strength of the biexciton resonance increases with uq. 



FIG. 2. The memory function S{uj) for attractive exciton-exciton interaction with a/l = 3.0 
and uq/Eq = —0.05 (solid line), uq/Eq = —0.4 (dotted line), and uq/Eq = —0.8 (dashed line). 
With increasing |no|, the asymmetric band of exciton-exciton scattering states transforms into a 
narrow antibound resonance. 



FIG. 3. Calculated TI-FWM signal for short pump duration, IoEq = 0.5, for a/l = 6.0, 
r/Eo = 0.0125, -yb/Eo = 0.005, Qq = 0, with uq/Eq = 0.05 (solid curve) and uq/Eq = 0.4 (dashed 
curve). The period of the biexciton oscillations is Ti, = 271 /E^ ~ 2tx /uq. 



FIG. 4. Calculated SR-FWM signal for uq/Eq = 0.4, and the rest of the parameters same as 
in Fig. p. The exciton peak strength oscillates with the biexciton period = 2TT/Eb ~ 27r/no. 



FIG. 5. Calculated TI-FWM signal at intermediate magnetic field, a/l = 3.0, for uq/Eq = 0.05 
(solid line), uq/Eq = 0.4 (dotted line), and uq/Eq = 0.8 (dashed line) with (a) short pump duration, 
toEo = 0.5, and (b) long pump duration, IqEq = 5.0. The rest of parameters same as in Fig. ^. 



FIG. 6. Calculated SR-FWM signal for uq = 0.4, pump tuned to exciton, = 0, and long 
pulse duration, IqEq = 5.0. The biexciton peak is suppressed while the exciton peak decays with 
characteristic time (2r)~^ for r > 0. The rest of parameters same as in Fig. |5|. 



FIG. 7. Calculated SR-FWM signal for the pump tuned to biexciton, 20,o = Ef,. For r > 0, 
biexciton peak decays with characteristic time (2r)~^ while exciton peak follows the pump. The 
rest of parameters same as in Fig. 0. 



FIG. 8. Comparison of TI-FWM signals calculated using the analytical solution ( p^pq ) (dotted 
curve and the APM (P) (dashed curve) with the exact numerical calculations (solid curve) for 
a/l = 1.5, uo/\Eo\ = 0.4 with (a) 76 = 0.0055^0 and (b) 7f, = 2r = 0.025£;o- The rest of parameters 
same as in Fig. H. 



FIG. 9. Calculated TI-FWM signal for repulsive interaction with G = —6.0 (solid line), 
G = -4.8 (dotted line), G = -3.6 (dashed line), G = -1.2 (long-dashed line), and G = -0.24 
(dot-dashed line), with (a) T/Eq = 0.0125 and (b) T/Eq = 0.05. For larger |G|, the decay is expo- 
nential with the interaction-induced characteristic time 7^^. For smaller the exponential decay 
is followed by the nonexponential tail for large |r|. The rest of the parameters are the same as in 
Fig. I 
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FIG. 10. (a) Comparison of TI-FWM signals calculated using analytical solution (p4t36|) (dot- 
ted curve) and the APM ( ^8[ ) (dashed curve) with the exact numerical calculations (solid curve) 
for a/l = 3.0, G = —1.2 with T /Eq and IqEq = 0.5. The exponential decay at initial time delays is 
followed by the power law decay for large |t|. (b) Real (upper curve) and imaginary (lower curve) 
parts of the memory function fitted by a Lorentzian (dotted line) . The Lorentzian width determines 
the exponential decay of the signal in the APM. The deviation of S{lo) from the Lorentzian at low 
frequencies leads to the power-law decay for large |r|. 



FIG. 11. Same as Fig. 10, but for weaker magnetic field, a/l = 1.5. The crossover to the 



nonexponential regime occurs at earlier time delays. 
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Fig. 6 



36 




37 




38 




-600 -400 -200 



Fig. 9 



39 




Fig. 10 



40 




41 



